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^ ! Abstract 
O 

The gluon propagator is calculated in quenched QCD for two different lattice 
sizes (16 3 x 48 and 32 3 x 64) at = 6.0. The volume dependence of the 
propagator in Landau gauge is studied. The smaller lattice is instrumental 
in revealing finite volume and anisotropic lattice artefacts. Methods for min- 
p imising these artefacts are developed and applied to the larger lattice data. 

New structure seen in the infrared region survives these conservative cuts to 
the lattice data. This structure serves to rule out a number of models that 
have appeared in the literature. A fit to a simple analytical form capturing 
the momentum dependence of the nonperturbative gluon propagator is also 
d ' reported. 
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I. INTRODUCTION 



The infrared behaviour of the nonperturbative gluon propagator has been the subject 
of extensive research and debate. Knowledge of this behaviour is generally regarded as 
being central to understanding the mechanism of confinement in quantum chromodynamics 
(QCD). Moreover, it is of great importance in various model-based approaches. For example, 
some studies based on Dyson-Schwinger equations have indicated that an infrared 

enhanced gluon propagator may be required for confinement; however, others @J5||| have 
pointed to the possibility of a dynamically generated gluon mass, or other forms leading to 
an infrared finite (or even vanishing) propagator (see [0] and references therein). 

Computing the gluon propagator directly on the lattice should provide an opportunity to 
resolve these contradictory claims. However, previous lattice studies ||[| have been unable 
to access the "deep" infrared region where the most interesting nonperturbative behaviour is 
expected. Significant finite volume effects introduced through zero momentum components 
prevent the study of the infrared behaviour of the propagator on a small lattice. 

The main aim of this study is to obtain a definite behaviour for the gluon propagator 
for momenta smaller than 1 GeV, where the interesting physics is expected to reside. In the 
following, we will report such results for momenta as small as 0.4 GeV. These momenta are 
small enough to reveal new structure in the gluon propagator and to place strong constraints 
on its infrared behaviour. 



II. THE GLUON PROPAGATOR ON THE LATTICE 

The gauge links U^x) may be expressed in terms of the continuum gluon fields as 

where V denotes path ordering. From this, the dimensionless lattice gluon field A^(x) may 
be obtained via 

Afc + A/2) = JL (y^x) - Ul(x)) - ^-Tr (y,{x) - E#(x)) , (2.2) 

which is accurate to 0(a 2 ). Denoting the discrete momenta available on a finite, periodic 
volume by q, the momentum space gluon field is 



A,(q)^J2e- i ^ + ^A^(x + il/2) 



A 4 

x 



—iq^a/2 r i 

" [(U,(q) - UH-q)) - -Tr (U,(q) - 



(2.3) 



where U^(q) = J2x e ^U^x) and A^(q) = t a A a ^(q). The available momentum values q are 
given by 



2 



& = -^» n M = 0,...,L M -l (2.4) 

where L M is the length of the box in the fi direction. The factor e _i ^ a / 2 , which comes from 
the gauge fields being defined on the links rather than the sites of the lattice, is crucial to 
removing (9(a)-errors and in particular obtaining the correct tensor structure for the gluon 
propagator |jID| . The gluon propagator D a ^ u {q) is defined by 

DfM = (A;(q)Al(-q)) . (2.5) 

We choose to study the gauge dependent propagator in the Landau gauge, which can be 
implemented numerically by maximising F [g] = J2i_i, x ^Tr (x) , where 

U°(x) =g(x)U,(x)g\x + fi). (2.6) 

In the continuum limit, this is related to the fact that fields satisfying the Landau gauge 
condition correspond to stationary points of F cont [g] =\\ A 9 || 2 = J dhcTr {A 9 ^) 2 {x) ||. 

The Landau gauge gluon propagator in the continuum has the form 

DtM = - ^)<*°W) , (2-7) 

The scalar function D(q 2 ) can then be extracted using 

^ 2 ) = ^E^E^)- (2-8) 

For a — > 0, the lattice propagator is related to the continuum one by a 2 D(q 2 ) = D cont (q 2 ) + 
0(a 2 ). 

A well-known lattice artefact is that the tree level propagator of a massless boson field 
does not reproduce the expected continuum result of 

£ ( °V) = ^, (2.9) 



but rather produces 



a 2 



a 2 D^\q) = — . (2.10) 

W E M (2sing>/2)2 1 } 



In the following, we are particularly interested in the quantity q 2 D(q 2 ), which is expected 
to approach 1 up to logarithmic corrections as q 2 —>■ oo. To ensure this result we will work 
with the momentum variable defined as§ 



x The momenta q and q are often defined the other way around in the lattice literature. However, 
we feel it is more instructive to define q as above, such that lattice results reproduce the continuum 
formula ( |2.7| ) and the tree level formula (2.9). 
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2 . %a 



(2-11) 



III. SIMULATION PARAMETERS 



We have analysed 75 configurations at (3 = 6.0, 



1.885 GeV based on the string tension in |TT 



on a 32 3 x 64 lattice. Using the value of 
this corresponds to a physical volume 



of (3.35 3 x 6.70)fm. For comparison, we have also studied an ensemble of 125 configurations 
on a smaller volume of 16 3 x 48, with the same lattice spacing. 

The gauge configurations were generated using a combination of seven over-relaxation 
and one Cabibbo-Marinari updates, with a separation between configurations of 1000 com- 
bined updates for the large lattice and 800 for the smaller lattice. Both lattices were fixed 
to Landau gauge using a Fourier accelerated steepest descent algorithm [O. An accuracy 



of vj^H^x |<VM 2 < 10 12 was achieved.! 



IV. FINITE SIZE EFFECTS AND ANISOTROPIC BEHAVIOUR 



A. Small Lattice Analysis 



We begin by considering the effect of the kinematic correction introduced through the 
change of variables in (|2.11|) . In the absence of this correction, data in the high momentum 
region are expected to display significant anisotropy when shown as a function of q. This is 
confirmed in Fig. [I], which shows the gluon propagator multiplied by q 2 a 2 and plotted as a 
function of qa. Here and in the following, a Z 3 averaging is performed on the data, where for 
example the momentum along (t,x,y,z) = (1,2,1,1) is averaged with (1,1,2,1) and (1,1,1,2). 

We expect the kinematic correction to reduce anisotropy, particularly in the large mo- 
mentum region. Fig. |2] shows the gluon propagator multiplied by the factor q 2 a 2 and plotted 
as a function of qa for momenta directed throughout the lattice. The anticipated reduction 
of anisotropy for qa > 1.5 is nicely displayed in this figure. A similar result was found in ||. 

Since the low momentum region holds the greatest nonperturbative interest, it is instruc- 
tive to stress that the low momentum points displayed in Figs, p] and ^| are insensitive to 
the whether one plots as a function of q or q. It is also useful to note that on any finite 
lattice, D(q 2 ) will be finite at q 2 a 2 = 0. Hence any lattice calculation must give q 2 a 2 D(q 2 ) 
vanishing at q 2 = 0. 

On the small lattice, we also see significant anisotropy in the data which have their origin 
in finite volume artefacts. Finite size effects are expected to be largest where the momentum 



2 The average link trace ^(KTrC/^} for the gauge-fixed configurations at /? = 6.0 was determined 
to 0.8609(7) for the small lattice and 0.8617(1) for the large lattice. 
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FIG. 1. The gluon propagator from our small lattice multiplied by q 2 a? plotted as a function of 
momenta qa. Values for each momentum direction are plotted separately. Only a Z3 averaging has 
been performed. Filled squares denote momenta directed along spatial axes, while filled triangles 
denote momenta directed along the time axis. Other momenta are indicated by open circles. 



component along one or more of the lattice axes is zero. Since the length of the lattice in 
the time direction is three times that of the spatial directions, we notice a clear difference 
between points which correspond to momenta directed along spatial axes from those purely 
in the time direction. These finite volume artefacts are clearly displayed at small momenta 
by the discrepancies between the filled squares (denoting momenta directed along spatial 
axes), and the filled triangles (denoting momenta directed along the time axis). 

Some residual anisotropy remains at moderate momenta around qa ~ 1.5 despite includ- 
ing the kinematical correction of ( 2.11 ). This anisotropy is clearly displayed in Fig. [|by the 
filled squares and triangles denoting momenta directed along lattice axes lying below the 
majority of points from off-axis momenta! for qa ~ 1.4. Since tree-level and one- loop 0(4) 
breaking effects should be removed by the kinematical correction, the remaining anisotropy 
appears to have its origin in quantum effects beyond one loop. 

In order to minimise lattice artefacts for large momentum components we select momen- 
tum vectors lying within a cylinder directed along the diagonal (t,x,y,z) = (1,1,1,1) of 



3 In Fig. I it appears that this anisotropy disappears as one goes to even larger momenta. However, 
this is not necessarily the case. Only momentum components up to q^a = 4 x 27r/16 have been 
selected in preparing Fig. §. This means that the largest momenta are not accessed through any 
single Cartesian direction. 
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FIG. 2. The gluon propagator from our small lattice multiplied by q 2 a? plotted as a function 
of momenta qa. Values for each momentum direction are plotted separately. Only a Z3 averaging 
has been performed. The symbols are as in Fig. |l]. Finite volume errors are expected to be 
largest for the purely time-like momenta, as the three shorter spatial directions have momentum 
components equal to zero and hence the effects of the nearby spatial boundaries are expected to 
be at a maximum. 



the lattice. This allows one to access the largest of momenta with the smallest of compo- 
nents. We found the selection of a cylinder with a radius of one spatial momentum unit 
(Aqa < 1x2tc/L s , where L s is the number of sites along a spatial axis) provides a reasonable 
number of points falling along a single curve for large momenta. The data surviving this 
cut are displayed in Fig. ||. 

However, this cut does not address the large finite volume errors surviving in Fig. ^. To 
remove these problematic momenta, we consider a further cut designed to remove momentum 
vectors which have one or more vanishing components. This is implemented by keeping only 
momentum directions that lie within a certain solid angle of the diagonal. We found that 
a cone of angle 20° measured from the diagonal at the origin was sufficient to provide a set 
of points lying along a smooth curve. Fig. |] displays these data. After these conservative 
cuts, there is relatively little structure left in the infrared region on our small lattice. Armed 
with this knowledge of how to obtain reliable lattice data, we now turn our attention to the 
gluon propagator data obtained from our larger lattice. 
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FIG. 3. The gluon propagator from our smaller lattice multiplied by q 2 a 2 . The points dis- 
played in this plot lie within a cylinder of radius Aqa < 1 x 2n/16 directed along the diagonal 
(t, x, y, z) = (1, 1, 1, 1) of the lattice. 
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FIG. 4. The gluon propagator from our smaller lattice multiplied by q 2 a 2 . The points displayed 
in this plot lie within a cylinder of radius Aqa < lx27r/16 directed along the diagonal of the lattice 
and within a cone of 20° measured from the diagonal at the origin. 
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FIG. 5. The gluon propagator from our larger lattice multiplied by q 2 a 2 plotted as a function 
of momenta qa. Values for each momentum direction are plotted separately. Only a Z3 averaging 
has been performed for the data shown in this figure. Plotting symbols are as in Fig. ||. Finite 
volume errors are greatly reduced compared to the results from the smaller lattice, as displayed by 
the overlap of points obtained from spatial and time-like momentum vectors. However, significant 
anisotropy is apparent for larger momenta. 



B. Large Lattice Analysis 



Fig. |^ displays the gluon propagator data for all momentum directions and values on 
the larger lattice. Again, only a Z 3 averaging has been performed. Examination of the 
infrared region indicates that finite volume artefacts are very small on the larger lattice. In 
particular, the agreement between purely spatial (filled squares) and time-like momentum 
vectors (filled triangles) at qa = 0.20 appears to indicate that finite size effects are relatively 
small on this lattice. Such an observation is consistent with topological studies of the QCD 
vacuum which provide some insight into the typical scale of QCD vacuum fluctuations. 

At large momenta, qa > 1.0, significant anisotropy is observed, similar to those displayed 
in Fig. [|. The fact that this anisotropy occurs at the same momentum values and with the 
same magnitude on both lattices confirms our previous argument that they result from finite 
lattice spacing errors as opposed to finite volume errors. Similar behaviour is expected in 
this region as the lattice action and lattice spacing are the same for our two lattices. To 
side-step this problem, we adopt the same cut as before. On this larger lattice, all momenta 
must lie within a cylinder of radius two spatial momentum units directed along the lattice 
diagonal. Fig. |6| displays the momenta surviving this cut. 

While the two points at qa = 0.20 indicate that finite volume errors are small on this 
lattice, we also consider the additional angular cut, requiring that all points must lie within 
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FIG. 6. The gluon propagator from our larger lattice multiplied by q 2 a 2 . The points displayed 
in this plot lie within a cylinder of radius Aqa < 2x27r/32 directed along the diagonal of the lattice. 
The first data point will be ignored in all fits, since it is not possible to assess the finite size effects 
for this point. The agreement between the next pair of data points indicates that finite size effects 
here are small. 
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FIG. 7. The gluon propagator from our larger lattice multiplied by q 2 a 2 . The points displayed 
in this plot lie within a cylinder of radius Aqa < 2x27r/32 directed along the diagonal of the lattice 
and within a cone of 20° measured from the diagonal at the origin. 



9 



a cone of opening angle 20° from the diagonal at the origin. The data surviving this cau- 
tiously conservative cut are illustrated in Fig. |?]. It is interesting to see that the turnover in 
q 2 a 2 D(q 2 ) in the infrared region survives even this extreme cut on our larger lattice. 



V. MODELLING THE PROPAGATOR 



The approach we have described here differs in several respects from that of |J. In 
particular, the momentum vectors surviving our cuts are very different. In || only a small 
number of low-lying spatial momentum values are used along with all available time mo- 
mentum values. In contrast, we have treated momenta in the spatial and time directions on 
an equal footing and selected momenta lying near the 4-dimensional diagonal, where lattice 
artefacts are expected to be minimal. In fact, nearly all the momenta used in || would be 
excluded by our cuts to the data. 

In addition, we include the kinematical correction of (2.11) . The authors of Ref . [0] do 



not include such a correction when performing their fits. Instead, their fits are constrained 
to the low-momentum region where such artefacts are hoped to be small. 

We have considered a number of models for the propagator which have been suggested 
in the literature, as well as some simple analytical forms which are intended to capture the 
essence of the nonperturbative gluon propagator. A detailed analysis of these models is 
currently in progress [|l"3|j . 

The data for the fit are those obtained on the large lattice with the cylindrical cut. To 
balance the sensitivity of the fit over the available range of qa, we average adjacent lattice 
momenta lying within Aqa < 0.005. In all fits the first point, at qa ~ 0.1, is omitted as it 
may be sensitive to the finite volume of the lattice. 

Our results so far indicate that the following analytical form 

D ^ q ^ = Z {(q 2 a 2 y+ a + (M 2 Y+ a + q 2 a 2 + M 2 ) ' ^ 

provides a satisfactory description of the data over a wide range of momenta. Our best fit 
to ( |5.1|) is illustrated in Figs. |8| and |[ This fit yields x 2 /dof = 3.5, a somewhat high value. 
However, if the first four points are omitted, a more acceptable value of 1.6 is found. 

We have studied the stability of the fits by varying the fitting range, and formula ( p . 1| ) 
turns out to give stable values for all the parameters over a wide region of fitting ranges. 
Our best estimates for the parameters, using all the available data, are 
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FIG. 8. The gluon propagator multiplied by q 2 , with nearby points averaged. The line illus- 
trates our best fit to the form defined in Eq. (|5.1| ). The fit is performed over all points shown, 
excluding the one at the lowest momentum value, which may be sensitive to the finite volume of 
the lattice. 



where the errors denote the uncertainties in the last digit (s) of the parameters. The first 
set of errors are statistical uncertainties in the fit to the full data set, and the second set of 
errors are due to fluctuations in the parameters as the fitting interval is varied. The estimate 
for M corresponds to a physical value for the parameter in the region of 700 MeV. 



VI. CONCLUSIONS 



We have performed an accurate, nonperturbative study of the infrared behaviour of the 
gluon propagator in the pure gauge theory on the lattice. We were able to isolate a set of 
data points for which systematic lattice errors seem to be negligible. Our data indicate an 
infrared behaviour less singular than 1/q 2 . This can be inferred from our plots by noticing 
the clear turnover in the behaviour of q 2 a 2 D(q) around q 2 = 1 GeV 2 . This is in agreement 
with suggestions from previous lattice results. In particular, there is evidence that the 
nonperturbative gluon propagator may be infrared finite. Our data appear to rule out 
models like those proposed in ||, which lead to an infrared enhanced propagator. 

Work in progress focuses on a detailed analysis of various analytic forms for the gluon 



propagator [ 13| . A stability analysis of the fit parameters is central to identifying the 
model best able to describe the gluon propagator in both the nonperturbative and the well 
known perturbative regimes. We are also exploring the possibility of extracting the value 
of D(q = 0) in the infinite volume limit, which one may be able to extrapolate from data 
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FIG. 9. The gluon propagator in physical units, with nearby points averaged. The line illus- 
trates our best fit to the form defined in Eq. fl5.l|) . The fit is the same as that shown in Fig. ||. 



The scale is taken from the value for the string tension quoted in [11 1. 



at different volumes. A study of the effects of Gribov copies |L4| and a complete analysis of 
the tensor structure of the gluon propagator are also issues under consideration. 

Finally, one promising line of research appears to be the study of the gluon propagator 
using improved lattice actions fl5j . These should yield a significant reduction of finite lattice 
spacing effects. Then, by performing simulations at larger values of the lattice spacing, one 
would be able to measure the propagator on larger physical volumes, thus gaining access to 
very low momentum values. For example, one may be able to use lattice spacings as large 
as 0.4 fm, so that a modest 16 4 lattice would have a physical size of 6.4 fm. 
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